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Abstract 

In this paper, the influence of aggregate size and volume fraction on shrink- 
age induced micro-cracking and permeability of concrete and mortar was in- 
vestigated. Nonlinear finite element analyses of model concrete and mor- 
tar specimens were performed. The aggregate diameter was varied between 
2 and 16 mm. Furthermore, a range of volume fractions between 0.1 and 0.5 
was studied. The nonlinear analyses were based on a 2D lattice approach 
in which aggregates were simplified as monosized cylindrical inclusions. The 
analysis results were interpreted by means of crack width and change of perme- 
ability. The results show that increasing aggregate diameter (at equal volume 
fraction) and decreasing volume fraction (at equal aggregate diameter) greatly 
increases permeability. 

1. Introduction 

Drying of cement based composites, such as concrete and mortar, induces 
cracking if shrinkage of the constituents is either internally or externally re- 
strained. For example, non- uniform drying leads to a moisture gradient, which 
results in non-uniform shrinkage of the specimen. Surface regions shrink faster 
than the inner bulk material, which results in surface cracking [lj]. Addition- 
ally, shrinkage might be restrained by aggregates within the composite 0. 
Aggregate-restrained shrinkage can lead to micro-cracking, which strongly in- 
fluences the transport properties of the material 0,0]- However, the evolution 
of micro-cracks and their dependence on the size and volume fraction of ag- 
gregates is not fully understood yet. In [9] it was observed that permeability 
increases with increasing aggregate size at a constant aggregate fraction. This 
result is surprising, since an increase of the aggregate size at a constant aggre- 
gate fraction is usually accompanied by a decrease in the volume of interfacial 
transition zones (ITZs), which are known to be more porous than the cement 
paste. One hypothesis is that an increase of aggregate diameter at constant 
aggregate fraction results in an increase of micro-crack width, which is closely 
related to permeability The objective of this work was to establish whether 
this size effect really occurs. This will undoubtedly enhance the understanding 



1 



of the link between micro-structure and macro-property, in particular the effect 
of microcracking on mass transport, which is a critical aspect for predicting 
durability and service-life. 

In the present work, shrinkage induced micro- cracking of concrete was analysed 
by means of the nonlinear finite element method. A lattice approach was 
used in combination with a damage-plasticity constitutive model, which was 
designed to result in mesh- independent responses |4|. This lattice approach is 
robust and computationally efficient. 

The present study is based on several simplifications. Shrinkage is represented 
by an eigenstrain, which was uniformly applied to the cement matrix only. This 
is representative of autogenous shrinkage, but does not fully represent tran- 
sient non-uniform shrinkage due to moisture gradients which may lead to more 
cracking near the surface than in the centre of a concrete element. Further- 
more, the only inclusions considered were aggregates, which were embedded 
in a uniform cement paste. Micro-cracking due to other inclusions, e.g. unhy- 
drated cement and calcium hydroxide crystals, was not considered. Aggregates 
were assumed to be separated from the cement paste by interfacial transition 
zones (ITZs), which were modelled to be weaker and more brittle than the 
cement matrix. Furthermore, the study was limited to two-dimensional plane 
stress analyses, in which aggregates are idealised as cylindrical inclusions of 
constant diameter. 

2. Modelling approach 

In the present work, shrinkage induced cracking was described by means of 
a lattice approach [3| combined with a damage-plasticity constitutive law jij. 
Nodes are placed randomly in the specimen constrained by a minimum distance 
d m , i.e. the smaller d m is, the smaller the average element length (Fig. 
Based on these randomly placed nodes, the spatial arrangement of lattice ele- 
ments is determined by a Voronoi tesselation. The cross-sections of the lattice 
elements, which connect the nodes (Voronoi sides), are the edges of the Voronoi 
polygons (Fig. [T^l) . Each node has three degrees of freedom, two translations 
and one rotation, shown in the local coordinate system in Figure [lb. The de- 
grees of freedom u e = {u\, </>i, u 2 , v 2 , </>2} T of two nodes of a lattice element 
are related to the displacement discontinuities u c = {u c , v c } T at the mid-point 
C of the cross-section. The constitutive model of the lattice elements is a 
combination of plasticity and damage mechanics [H, which relates stresses to 
strains. The strains e = (e n ,e s ) T are determined from the displacement jump 
u c = (un,u s ) T at mid-point C as 




Figure 1: Discretisation: (a) Lattice elements (solid lines) and cross-sections 
(dashed lines) obtained from the Delaunay triangulation and dual Voronoi 
tessellation, respectively, (b) Lattice element with degrees of freedom in the 
local co-ordinate system. 

These strains are related to the nominal stress cr = (cr n , a s ) T as 

a = (1 - u/)D e (e - e p - e s ) = (1 - u) & (2) 

where u> is the damage variable, D e is the elastic stiffness, e p = (s pn ,e pB ) T is 
the plastic strain, s s = (e s ,0) is the shrinkage eigenstrain and & = (a n ,a s ) 
is the effective stress. 

The plasticity part of the model is based on the effective stress a and consists of 
an elliptic yield surface, an associated flow rule, an evolution law for the hard- 
ening parameter and loading and unloading conditions. A detailed description 
of the components of the plasticity model is presented in jij. The initial yield 
surface of the plasticity model is determined by the tensile strength f t , the 
shear strength s/ t and the compressive strength cf t . The evolution of the 
yield surface during hardening is controlled by the model parameter fi, which 
is defined as the ratio of permanent and reversible inelastic displacements. The 
damage part is formulated so that linear stress inelastic displacement laws for 
pure tension and compression are obtained, which are characterised by the 
fracture energies Gf t and Gf c , respectively. 

The equivalent crack opening w c is defined by the equivalent inelastic displace- 
ment, which for the present damage-plasticity model is defined as 

w c = || w c || (3) 
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Figure 2: Analysis results at a shrinkage strain of e s = 0.5 %: (a) Deformations 
(magnified by a factor of 50) for <fi = 16 mm. (b) Crack permeability k c (Eq. [6]) 
versus aggregate diameter <fi for volume fractions of p = 0.1, 0.3 and 0.5. 

where 

w c = h(e p + u (e - e p )) (4) 

The inelastic displacement vector w c is composed of a permanent and reversible 
part, defined as he p and hu> (e — e p ), respectively. 

3. Nonlinear finite element analysis of shrinkage induced micro- 
cracking 

Shrinkage induced micro- cracking was analysed by means of the nonlinear 
lattice approach described above. The elements representing the cement paste 
were subjected to an incrementally applied uniform shrinkage strain up to 
e s = 0.5 % (Eq. |2J). This value was chosen for the simulation to represent 
a relatively severe shrinkage of the neat cement paste on first-drying to low 
humidities. The influence of aggregate volume fraction and aggregate diameter 
was studied. Aggregate volume fractions p = 0.5,0.3 and 0.1 were modelled. 
Furthermore, four different aggregate diameters ((f) = 16, 8, 4, and 2 mm) were 
used. The length L of the specimen 




Thus, the smaller the aggregate diameter and the greater the volume ratio 
p, the smaller is the specimen length L (Fig. [2]). The smallest separation dis- 
tance between aggregate particles, i.e. minimum width of the cement paste 
decreases with increasing aggregate volume fraction at constant aggregate size, 
as would be expected in the case of real mortars with increasing aggregate frac- 
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tion. At constant aggregate fraction, the separation between aggregate parti- 
cles increases with increasing aggregate size, similar to the case of a mortar 
compared to concrete at the same aggregate content. Aggregates were mod- 
elled elastically. Lattice elements crossing the boundary between aggregates 
and cement paste represent the average response of the interfacial transition 
zones and the two adjacent material phases. In all the analyses, the length 
of the lattice elements is significantly greater than the width of the highly 
non-uniform interfacial transition zones, which is usually in the range of 10s 
of fim [3]. Therefore, the stiffness of these lattice elements was approximated 
as an average of the Young's modulus of aggregate and mortar. The strength 
of these lattice elements is determined by the strength of the ITZ, which is 
the weakest link. Here, the strength and fracture energy ratio of cement paste 
and ITZ was chosen as 2. This ratio is an approximate value for samples with 
a relatively weak ITZ. For an aggregate diameter of = 16 mm and a volume 
fraction of p = 0.3, the deformed mesh is shown in Fig. [2^. The deformations 
are localised in a few cracks, which connect the aggregates in a regular square 
pattern. Qualitatively very similar crack patterns were obtained for the other 
aggregate diameters and volume fractions. This pattern of cracking seems 
consistent to that observed by Hsu |6] in '2D' model samples made of sand- 
stone discs arranged in a square grid and filled with paste that is subsequently 
subjected to drying shrinkage. Depending on the separation between the ag- 
gregates, cracks were seen to occur at the interface (i.e. bond cracks), near the 
shortest distance and diagonally at the largest distance between aggregates. 

The average crack width increases with increasing aggregate size at constant 
aggregate volume fraction. At constant aggregate size, the average crack width 
increases with decreasing aggregate volume fraction. Crack width is closely 
related to transport properties of concrete, in particular in the case of flow 
under a pressure gradient. Assuming that the paste matrix is dense so that 
flow occurs predominantly through the cracks, the permeability k c in the out- 
of-plane direction due to cracking can be described by the cube of the crack 
width as 



where £ is a material constant [8(]. The influence of aggregate diameter and 
volume fraction on k c is shown in Fig. [2b on a log- log scale for £ = 1. Note that 
what is of interest here is not the actual value of the estimated permeability, 
but the change in permeability caused by varying either the aggregate volume 
fraction or particle size. At a constant aggregate volume fraction, increasing 
the aggregate diameter from 2 to 16 mm caused approximately a 2.5 orders of 
magnitude increase in permeability. At constant aggregate diameter, increas- 
ing the aggregate volume fraction from 0.1 to 0.5 produced about 1 order of 
magnitude decrease in permeability. The aggregate diameter influences the 
permeability strongly, since the crack width, which increases with increasing 




(6) 



8=1 
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aggregate diameter, enters Eq. (jSj) in its cube. 
4. Conclusions 

In the present work the influence of aggregate size and volume fraction on 
shrinkage induced micro-cracking was studied numerically by means of the 
nonlinear finite element method. The work has shown that the average crack 
width increases with increasing aggregate diameter and decreasing volume 
fraction. The permeability, which is related to the cube of the crack width, 
increases with increasing aggregate diameter and decreasing volume fraction 
Thus, the aggregate diameter and volume fraction influence the formation of 
microcracks induced by aggregate restrained shrinkage, which is in agreement 
with experimental observations in [9J. However, the representation of aggre- 
gates as circular inclusions and the use of one aggregate size oversimplifies the 
micro-structure of concrete, which does not allow a direct comparison with 
experimental results. In future work, the present modelling approach will be 
extended to 3D and the use of realistic aggregate grading curves will be in- 
cluded by applying multi-scale analysis approaches. 
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